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ABSTRACT 


A  constitutive  model  for  ceramics  by  Johnson  and  Holmquist  has  been  incorporated 
and  validated  in  the  finite  element  computer  program  LS-DYNA2D.  Using  experimental 
data  acquired  at  DREV,  parameters  for  high  density  ceramics  (around  98%  theoretical 
maximum  density)  have  been  established.  The  effect  of  varying  some  of  the  parameters 
on  the  calculated  depths  of  penetration  into  semi-infinite  targets  and  the  ballistic  limits  of 
armour  has  been  investigated. 


RESUME 


Le  modele  de  ceramique  de  Johnson  et  Holmquist  a  ete  utilise  et  verifie  dans  le 
programme  d’elements  finis  LS-DYNA2D.  Utilisant  des  donnees  experimentales  du  CRDV, 
les  parametres  pour  des  ceramiques  de  haute  densite  (environ  98  %  de  densite  maximale) 
ont  ete  etablis.  L’effet  des  variations  de  ces  parametres  sur  les  profondeurs  de  penetration 
dans  des  cibles  semi-infinies  et  sur  les  limites  balistiques  a  ete  etudie. 
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EXECUTIVE  SUMMARY 


Current  emphasis  in  numerical  simulations  is  on  accurate  modelling  of  material 
behaviour.  Modelling  is  cheaper,  safer  and  more  flexible  than  doing  many  series  of  experi¬ 
ments.  However,  the  computational  description  of  a  material  must  be  as  accurate  as  possible 
before  simulations  of  it  are  considered  reliable.  Ceramics  and  other  brittle  materials  are 
difficult  to  model  and  characterize  because  damage  and  fracture  mechanisms  are  not  easily 
definable. 

This  memorandum  describes  one  model,  the  Johnson-Holmquist  model,  which  was 
incorporated  into  a  finite  element  hydrodynamic  computer  program  called  LS-DYNA2D. 
The  incorporation  was  verified  against  previously  reported  work  and  was  found  to  be  cor¬ 
rectly  incorporated.  Simulations  of  impacts  of  projectiles  (SLAP,  APM-2)  onto  high  density 
ceramic  tiles  covering  aluminum  were  performed.  Various  thicknesses  of  ceramic  were  used, 
corresponding  to  those  used  in  depth  of  penetration  and  ballistic  limit  tests  at  DREV.  By 
adjusting  the  values  of  key  parameters  in  the  model,  these  simulations  were  fitted  to  agree 
with  experimental  data. 

With  the  model  and  these  fitted  values,  DREV  is  in  a  better  position  to  simulate 
ceramic  behaviour  at  high  strain  rates.  Such  simulations  will  be  of  special  importance  in 
the  ballistic  protection  area,  for  both  personnel  and  vehicles  (LAV,  tanks)  where  full  scale 
tests  are  very  expensive. 
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LIST  OF  SYMBOLS 
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intact  material  strength  parameter 
damaged  material  strength  parameter 
strain  rate  influence  parameter 
damage  parameters 
failure  strain 
shear  modulus 

elastic  constants  for  Johnson-Holmquist  model 
damaged  material  strength  parameter 
intact  material  strength  parameter 
Pressure  (may  be  subscripted) 
tensile  strength 

fraction  of  damage  energy  to  convert  to  internal  pressure 

strain 

p/po~  1 

density  (may  be  subscripted) 
stress 

superscript  indicating  a  normalised  quantity 

depth  of  penetration 

Hugoniot  elastic  limit 

light  armoured  vehicle 

Livermore  Software  Technology  Corporation 

rolled  homogenous  armour 

saboted  light  armour  piercing 

theoretical  maximum  density 


The  above  list  is  not  exhaustive  but  does  include  most  of  the  symbols  used  in  this 


document. 


UNCLASSIFIED 

1 


1.0  INTRODUCTION 

One  of  the  current  activities  in  numerical  simulations  is  the  modelling  of  material 
behaviour  under  high  velocity  impact  since  it  is  cheaper,  safer  and  more  flexible  than  doing 
many  series  of  experiments.  However,  the  computational  description  of  a  material  must 
be  as  accurate  as  possible  before  simulations  of  it  can  be  considered  reliable.  One  type 
of  material  which  is  of  current  interest,  especially  for  armour  applications,  is  ceramics. 
A  ceramic  material  has  a  complicated  response  to  high  strain  rates,  with  damage  and 
fracture  being  very  important.  One  semi-empirical  model  to  describe  the  behaviour  of  brittle 
materials  such  as  ceramics  under  high  dynamic  loads  is  called  the  Johnson-Holmquist  model 
(Ref.  1,2).  In  this  memorandum,  this  model  will  be  described  and  its  implementation  and 
validation  in  the  finite  element  hydrodynamic  computer  program  called  DYNA2D  (Ref.  3) 
will  be  given. 

This  work  was  performed  at  DREV  between  October  1996  and  April  1998  under 
Work  Unit  2bbl5,  Numerical  Modelling  of  Ballistic  Events. 

2.0  MODEL  DESCRIPTION  AND  INCORPORATION  TESTS 

The  Johnson-Holmquist  model  was  originally  formulated  (Ref.l)  to  consider  the 
behaviour  of  brittle  materials  and  then  further  refined  (Ref.  2)  to  be  more  flexible  and 
less  parameter  sensitive.  The  constitutive  behaviour  described  by  this  model  is  as  follows. 
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There  is  an  initial  elastic  regime  which  continues  until  yield  occurs  (the  plastic  strain  starts 
at  this  point).  Damage  starts  to  accumulate  and  the  material  weakens.  The  material  then 
behaves  along  the  weakened  material  curve.  The  rate  at  which  the  material  moves  from 
undamaged  to  damaged  and  the  strength  of  the  the  weakened  material  are  all  variables  the 
user  must  establish  but  unfortunately,  they  are  not  directly  measurable.  Hence,  simulations 
with  various  parameters  are  necessary  to  establish  a  self-consistent  set  useful  for  describing 
the  material’s  behaviour. 

The  model  will  now  be  described  in  some  detail.  In  the  following  equations,  the 
stresses  and  pressures  are  normalized  as  follows.  Normalized  quantities  are  indicated  by  an 
asterix  superscript.  Given  the  Hugoniot  elastic  limit  (HEL)  and  the  shear  modulus  (G), 
/ 1(=  pi p0  —  1)  at  the  HEL,  phel,  is  found  by  solving 

HEL  =  kipuEL  +  k2  Peel  +  ^3  Phel  +  — 

6  1  + PHEL 

where  ki,k2,  k$  are  model  elastic  constants.  Next,  the  pressure  at  the  HEL,  Phel,  is  found 
via 

Phel  —  kiPHEL  +  k2pjjEL  +  ^3  P3hel 
and  the  stress  at  the  HEL,  &hel,  from 

&hel  =  HEL  -  Phel)- 
Then,  the  normalized  pressure,  P* ,  is  given  by 


P*  =  P/ Phel 
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where  P  is  the  real  pressure.  Any  normalised  stress,  a*,  is  given  by 

a*  -  a/cruEL 

where  a  is  the  real  stress.  The  description  of  the  model  can  now  begin.  In  a  brittle  material 
subjected  to  damage,  the  equivalent  stress  is 

a*  =  of  ~  D(<r*  -  a}) 

where  the  strength  envelope  for  the  intact  material  is 

a*  =  a(P*  +  T*)n(l  +  cine), 
the  strength  of  the  damaged  material  is 

c}  =  b(P*)m(l  +  c\ne), 

and  the  damage  parameter  is 

D  =  £  Aep/j, 

which  depends  upon  the  plastic  strain  to  fracture: 

=  d\(P*  +  T*)di 

and  the  incremental  plastic  strain  in  a  given  computational  cycle,  .  In  the  above, 
a,n,c,b,m,di,d2  are  all  material  constants,  T*  is  the  normalised  tensile  strength  and  e  is 
the  normalised  strain  rate.  As  for  the  pressure,  it  is  given  by 

P  =  kifj,  +  fc2//2  +  k3fi3  +  AP 
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(in  compression)  or  by 


P  =  kifi 


(in  tension)  where  A P  is  a  pressure  increment  due  to  the  conversion  of  energy  loss  (due  to 
damage)  into  internal  energy  (pressure).  Note  that,  in  this  model,  any  frictional  energy  loss 
is  considered  negligible. 


The  Johnson-Holmquist  model  was  incorporated  into  LS-DYNA2D  (Ref.  3)  by  the 
company  Livermore  Software  Technology  Corporation  (LSTC)  and  corrected  at  DREV.  A 
test  case  (described  in  Ref.  2)  was  simulated.  It  consists  of  a  controlled  compression  at  a 
constant  rate  followed  by  a  constrained  release  at  the  same  rate  (opposite  direction).  The 
resulting  stress-pressure  curve  is  shown  in  Fig.  1  as  the  solid  curve.  Various  features  of 
the  model  are  visible  in  this  figure.  There  is  an  initial  linear,  elastic  behaviour  until  the 
a*  curve  (solid  circle  symbols)  is  reached.  Damage  starts  at  this  point  and  continues  to 
increase  until  it  reaches  the  aj  curve  (solid  square  symbols).  As  the  applied  force  direction 
is  reversed,  elastic  release  occurs.  As  the  system  goes  through  its  initial  position,  the  elastic 
force  starts  to  rise  again  (due  to  the  internal  pressure  resulting  from  damage)  and  reaches 
the  limiting  aj  curve  which  it  follows  again.  This  behaviour  matches  the  results  given  in 
Ref.  2  and  were  reproduced  with  only  one  parameter  change:  our  dx  was  0.00815  instead 
of  original  0.005.  This  gives  one  confidence  that  the  model  is  incorporated  correctly.  The 
plastic  strain  as  calculated  in  LS-DYNA2D  may  be  different  to  that  of  Ref.  2  in  which 
the  finite  element  program  EPIC  was  used.  EPIC  uses  triangular  finite  elements  and  must 
therefore  handle  the  strains  differently.  Thus,  the  damage  may  accumulate  differently  as 
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TABLE  I 


Johnson-Holmquist  parameters  for  an  around  98%-TMD  ceramic 


Parameter 

Value 

P 

G 

1.5196  MBar 

*1 

2.30974  MBar 

-1.60027  MBar 

k3 

23.7376  MBar 

to 

HEL 

0.0657  MBar 

T 

0.00262  MBar 

a 

0.88 

n 

0.64 

c 

0.07 

b 

0.28 

m 

0.6 

® max 

1.0 

0 

1.0 

d\ 

0.01 

d,2 

0.7 

is 

1.5 

well.  If  this  is  indeed  the  case  and  di  being  a  damage  parameter,  it  follows  that  di  would 
have  to  be  “code-dependent”  and  thus  different  in  the  two  cases. 


For  all  subsequent  simulations,  a  different  set  of  parameters  was  used  which  de¬ 
scribed  better  the  around  98%  theoretical  maximum  density  (TMD)  ceramic  material  that 
was  used  experimentally.  These  parameters  are  given  in  Table  I  and  were  originally  found 
in  the  code  EPIC.  In  the  parameter  variation  studies  which  follow,  these  parameters  were 
unchanged  except  where  specifically  noted.  In  general,  only  b  and  fs  were  changed. 


Stress  (Pa) 
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Pressure  (Pa) 


FIGURE  1  -  Ceramic  material  response  to  compression  and  release,  after  Ref.  2 
with  di=0. 00815.  The  solid  line  is  generated  by  the  application  and 
release  of  a  load. 
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3.0  VALUES  FOR  JOHNSON-HOLMQUIST  MODEL  PARAMETERS 

An  isolated  model  is  not  worth  very  much  in  the  absence  of  values  for  the  parameters 
within  it.  To  determine  these  values,  comparisons  of  simulations  with  experimental  results 
are  used,  the  values  being  varied  until  good  agreement  between  simulation  and  experiment 
is  reached.  At  DREV,  two  types  of  experimental  results  are  available  for  comparison:  depth 
of  penetration  tests  and  ballistic  limit  tests.  The  methods  used  and  results  obtained  from 
these  tests  will  be  described  briefly  in  the  following  subsections. 

A  depth  of  penetration  (DOP)  test  is  one  in  which  a  sample  of  a  material  is  af¬ 
fixed  to  a  large  block  of  backing  material  such  as  aluminum  or  rolled  homogeneous  armour 
(RHA)  and  a  projectile  strikes  the  top  surface  of  the  combined  face  plate-backing.  The 
resulting  hole  or  depression  in  the  backing  material  is  then  measured,  usually  as  a  func¬ 
tion  of  sample  material  thickness.  For  comparison  purposes,  the  experimental  depths  of 
penetration  obtained  at  DREV  were  5.63  cm  for  the  0.82  cm  thick  ceramic,  1.68  cm  for 
the  1.27  cm  thick  ceramic  and  0.83  cm  for  the  2.54  cm  thick  ceramic.  By  varying  certain 
material  parameters  in  the  simulations,  the  calculated  depth  of  the  hole  can  be  adjusted 
to  match  the  experimental  results.  In  these  simulations,  the  projectile  was  a  saboted  light 
armour  piercing  (SLAP)  tungsten  core  and  the  backing  material  was  6061-T6  aluminum, 
thick  enough  (usually  7  cm)  so  that  rear  surface  or  boundary  effects  did  not  affect  the  results 
(“semi-infinite”).  The  overlying  ceramic  facing  was  one  of  three  possible  thicknesses:  0.82 
cm,  1.27  cm  or  2.54  cm.  The  meshing  was  uniform  in  the  vertical  directions  within  each 
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material  and  variable  horizontally,  with  finer  meshing  at  the  center  and  coarser  meshing 
in  the  outer  regions  away  from  the  projectile’s  point  of  impact.  Finer  meshing  is  used  to 
resolve  better  the  interaction  between  the  small  tip  of  the  projectile  and  the  much  larger 
target. 

In  a  ballistic  limit  test,  a  plate  of  potential  armour  material  is  affixed  to  a  backing 
plate  of  finite  thickness.  The  velocity  of  the  impacting  projectile  is  adjusted  so  that  the 
projectile  is  just  stopped  (or  barely  perforates  the  target).  This  velocity  is  called  the 
ballistic  limit  and  is  closely  related  to  the  experimental  V50,  the  velocity  at  which  50%  of 
the  projectiles  will  perforate  the  target.  Experimentally,  random  variations  from  shot  to 
shot  (e.g.  in  yaw,  velocity,  alignment,  material  properties,  etc.)  cause  variable  results,  hence 
a  probabilistic  velocity  is  usually  reported.  In  the  present  simulations,  the  projectile  was 
an  APM-2  0.50  caliber  Mo-Mn  steel  core  and  the  backing  material  was  2024-T3  aluminum 
0.7225  cm  thick  while  the  ceramic  face  plate  was  0.8  cm  thick.  Uniform  meshing  was  used 
vertically  and  variable  meshing  horizontally,  with  finer  meshing  near  the  projectile’s  impact 
point  and  coarser  meshing  further  out  radially. 

There  are  many  parameters  in  the  Johnson-Holmquist  model  (see  Table  I)  and  it 
is  not  wise  to  try  to  fit  them  all.  That  would  take  too  much  time  and  tends  to  detract 
from  the  underlying  question  of  which  parameters  are  important  and  which  are  not.  The 
following  philosophy  was  adopted  in  this  document.  As  far  as  the  material  properties  of 
the  ceramic  are  concerned,  it  is  assumed  that  the  properties  of  the  undamaged  ceramic  are 
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known  best  and  hence,  can  be  used  without  modification  directly  as  found  in  the  literature 
(e.g.  in  the  code  EPIC)  or  can  be  determined  experimentally.  However,  for  the  damaged 
material,  freedom  is  allowed  and  the  different  values  can  be  tried.  Upon  examination  of  the 
equations  describing  the  ceramic  model,  various  critical  parameters  are  noted.  For  damage 
production,  the  critical  parameter  is  di.  For  the  strength  of  the  damaged  materials,  one 
parameter,  namely  b,  is  of  great  importance.  These  are  the  main  parameters  which  will  be 
examined  in  the  following  sections. 

3.1  Failure  Criteria 

In  a  penetration  problem,  since  there  are  no  available  fracture  models,  there  must 
be  some  mechanism  for  removing  material  for  penetration  to  occur.  Otherwise,  the  material 
will  simply  be  stretched,  sometimes  beyond  physical  limits  and  worse,  penetration  will  never, 
occur.  In  this  section,  two  items  will  be  addressed.  First,  the  most  common  erosion  criterion, 
a  critical  failure  strain,  will  be  discussed  with  a  view  to  finding  a  rational  approach  to 
selecting  an  appropriate  failure  strain.  This  is  the  criterion  used  for  most  of  the  subsequent 
simulations  in  this  document.  Second,  a  more  physical  criterion  is  strongly  desirable  and 
thus  other  possible  failure  criteria  are  outlined. 

One  common  erosion  criterion  used  in  hydrocodes  is  based  upon  a  failure  strain 
parameter  (f4):  an  element  is  eroded  if  the  plastic  strain  in  it  exceeds  a  user-specified  limit. 
However,  a  rational  approach  to  determining  a  failure  strain  cannot  be  done  independently 
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of  other  parameters  as  the  strains  generated  within  a  material  depend  upon  them.  One 
critical  parameter  for  the  J-H  model  is  b,  which  determines  the  strength  of  the  damaged 
material  and  it  should  be  determined  first.  One  starts  with  preliminary  values  for  all  the 
parameters,  taken,  for  example,  from  Table  I  and  uses  a  large  failure  strain  (for  example, 
10),  far  above  any  plastic  strain  generated  in  the  material  (and  hence  erosion  is  suppressed). 
A  simulation  of  a  DOP  test  is  done  and  the  magnitudes  of  the  resulting  strains  are  noted. 
A  new,  lower  failure  strain,  one  that  ensures  erosion,  is  then  chosen  and  the  simulation 
is  rerun.  The  DOP  is  noted  and  b  is  varied  to  get  agreement  between  experiment  and 
simulation.  At  this  point,  one  should  have  a  ‘good’  set  of  parameters  which  characterize 
the  material  according  to  the  constitutive  model  described  earlier.  The  ballistic  limit  tests 
are  then  simulated  with  this  set  of  parameters.  The  failure  strain  is  adjusted  so  that  with 
the  experimental  ballistic  limit  velocity,  a  projectile  will  barely  perforate  the  target.  The 
transition  to  penetration  is  quite  sensitive  to  the  value  of  fs.  In  Fig.  2,  the  residual  velocity 
(either  perforation  or  rebound)  is  shown  as  a  function  of  fs  for  two  different  b’s.  As  can  be 
seen,  there  is  a  steep  slope  around  a  critical  f5.  With  the  new  critical  failure  strain  (usually 
lower),  the  DOP  is  reconfirmed.  In  effect,  what  one  is  doing  is  iterating  to  fit  both  the  DOP 
and  ballistic  limit  tests  by  adjusting  two  parameters,  b  and  fs.  The  net  result  must  be  a 
reasonable  agreement  between  experiment  and  simulations  in  both  cases. 

Other  erosion  criteria  can  be  envisaged.  For  example,  one  is  a  damaged  strength 
model  in  which  a  fully  damaged  material  in  tension  beyond  its  (intact)  tensile  strength 
is  eroded.  This  is  relatively  easy  to  implement  and  was  done  so  by  eroding  any  element 


Residual  velocity  (m/s) 
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Failure  strain 


FIGURE  2  -  Variation  of  residual  velocity  with  failure  strain  for  an  impact  velocity 
of  400  m/s  onto  a  ceramic- aluminum  sandwich 


UNCLASSIFIED 

12 


with  P*  +  T*  less  than  0  and  for  which  the  damage  parameter  was  1  (fully  damaged). 
The  results  indicate  that,  without  changing  other  parameters  of  the  model,  the  material 
tends  to  erode  too  quickly  and  hence  DOP  and  ballistic  limit  results  vary  significantly  from 
experiments.  As  in  the  failure  strain  case,  varying  the  strength  of  the  damaged  material 
by  changing  b  results  in  better  agreement.  This  was  done  and  the  resulting  DOP  versus 
b  for  various  thicknesses  of  ceramic  are  shown  in  Fig.  3.  As  the  strength  of  the  damaged 
material  increases  (larger  b),  the  DOP  decreases.  The  experimental  results  cannot  be 
reproduced  exactly  but  using  a  b  around  0.75  yields  acceptable  results.  For  comparison, 
the  undamaged  material  has  a  strength  parameter  (a)  of  0.88  and  thus,  it  is  seen  that  to 
reproduce  the  experiments,  with  this  erosion  criterion,  the  material  cannot  lose  a  lot  of 
strength  when  damaged,  something  which  is  not  the  case  in  reality.  For  this  reason,  further 
investigation  of  this  criterion  was  not  pursued. 

3.2  Depth  of  Penetration 

For  these  simulations,  the  failure  strain  erosion  criterion  was  used.  Various  failure 
strains  were  tried,  keeping  the  damaged  material  strength  parameter  constant  (b=0.43).  In 
Fig.  4,  the  DOP  is  shown  as  a  function  of  failure  strain  and  for  3  different  thicknesses  of 
ceramic.  As  expected,  as  the  failure  strain  limit  is  raised,  the  DOP  decreases.  This  occurs 
because  the  ceramic  material  resists  the  penetration  longer  at  higher  failure  strains  and 
hence  slows  the  projectile  more  and  a  slower  projectile  does  not  penetrate  as  far  into  the 
infinite  backing  material.  A  failure  strain  of  1.5  yields  results  close  to  the  experimental 
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FIGURE  3  -  Variation  of  depth  of  penetration  with  strength  of  damaged  material 
using  a  pressure-based  erosion  criterion 
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FIGURE  4  -  Depth  of  penetration  into  aluminum  behind  ceramic  plates  of  various 
thicknesses  as  a  function  of  assumed  failure  strain  and  for  b=0.43 

values  but  it  is  impossible  to  find  a  fs  for  which  all  are  reproduced.  Near  fs  of  0.7±0.1,  a 
sharp  rise  in  the  DOP  is  noted.  This  value  of  fs  is  suggestive  as  it  is  near  the  f5  needed  to 
reproduce  the  ballistic  limit  results  as  discussed  in  the  previous  section. 


Using  the  failure  strain  determined  above  (1.5),  variations  of  the  damaged  mate¬ 
rial  strength  parameter  (b)  were  tried.  These  results  are  shown  in  Fig.  5  where  they  are 
compared  with  experimental  data  (dashed  line).  The  experimental  results  are  bracketed 
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by  these  variations  (less  than  ±  10%).  However,  none  of  the  three  b’s  reproduced  the 
experimental  results  exactly.  The  b  of  0.43  is  a  reasonable  compromise. 

The  final  set  of  simulations  was  done  with  fs=1.5  and  b=0.43,  but  the  damage 
parameter  dj  was  varied.  This  parameter  controls  the  rate  at  which  damage  accumulates. 
If  di  is  zero,  the  damage  is  ‘instantaneous’.  These  were  done  for  the  three  different  ce¬ 
ramic  thicknesses  and  the  results  are  plotted  in  Fig.  6.  Once  again,  substantial  variations 
were  observed  but  the  experimental  results  could  not  be  reproduced.  The  original  damage 
parameter  of  0.01  gives  reasonable  results  but  using  0.0125  would  give  better  agreement. 

3.3  Ballistic  Limit  Simulations 


For  these  simulations,  the  failure  strain  erosion  criterion  was  used.  Ballistic  lim¬ 
its  are  determined  in  simulations  by  plotting  the  residual  velocity  (velocity  with  which  a 
projectile  exits  the  target)  versus  the  impact  velocity.  This  was  done  for  three  b-f5  combi¬ 
nations  and  the  curves  are  given  in  Fig.  7  along  with  the  experimental  results.  For  b=0.43 
and  fs=1.5,  the  combination  which  gave  the  compromise  results  in  the  DOP  penetration 
case,  there  is  a  significant  difference  between  the  experimental  and  numerical  results.  By 
using  fs=0.6  (determined  to  be  the  best  for  the  ballistic  limit  velocity),  excellent  agreement 
is  obtained  at  all  impact  velocities.  An  intermediate  result  is  obtained  with  fs=1.5  and 
b=0.40,  a  combination  which  gives  an  upper  limit  to  the  results  in  the  DOP  case.  The  need 
for  two  quite  different  fs’s  to  reproduce  the  results  from  two  different  experiments  suggests 


Thickness  of  ceramic  (cm) 


FIGURE  5 


Depth  of  penetration  into  aluminum  versus  thickness  of  ceramic  facing 
for  various  damaged  material  strengths.  A  failure  strain  of  1.5  was 
used.  The  dashed  line  is  the  experimental  data. 
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Damage  parameter  (dl) 


FIGURE  6  -  Depth  of  penetration  into  aluminum  versus  the  damage  parameter,  di 
for  three  different  ceramic  thicknesses 
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Impact  velocity  (m/s) 

FIGURE  7  -  Comparison  of  calculated  ballistic  limits  with  experimental  data 

that  the  failure  strain  criterion  is  not  the  best  criterion  for  determining  erosion.  Further 
study  is  needed  here. 


4.0  CONCLUSIONS  AND  RECOMMENDATIONS 


Two  sets  of  parameters  were  found  which  can  be  used  to  model  around  98%  TMD 
ceramics  in  depth  of  penetration  tests  (Table  I  with  fs=1.5  and  b=0.43)  and  ballistic  limit 
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tests  (Table  I  with  fs=0.6  and  b=0.43).  More  work  is  needed  to  derive  a  better  erosion 
criterion  and  to  be  able  to  use  one  set  of  parameters  (perhaps  Table  I  with  f5=0.8  and 
b=0.43)  for  all  simulations  involving  this  ceramic  material. 


UNCLASSIFIED 

20 


5.0  REFERENCES 

1.  Johnson,  G.R.,  Holmquist,  T.J.,  “A  computational  constitutive  model  for  brittle  mate¬ 
rials  subjected  to  large  strains,  high  strain  rates  and  high  pressures”,  in  Shock- wave  and 
High-strain  Rate  Phenomena  in  Materials,  pp. 1075-1081,  editors  Meyers,  M.A.,  Murr, 
L.E.,  Staudhammer,  K.P.,  Marcel  Dekker  Inc,  New  York,  1992 

2.  Johnson,  Gordon  R.,  Holmquist,  Tim.  J,  “An  improved  computational  constitutive 
model  for  brittle  materials”,  in  High-Pressure  Science  and  Technology  -  1993,  American 
Institute  of  Physics  Conference  Proceedings  No.  309,  pp. 981-984,  editors  Schwartz,  S.C., 
Shaner,  J.W.,  Samara,  G.A.,  Ross,  M.,  AIP  Press,  New  York,  1994 

3.  Hallquist,  John  0.,  “LS-DYNA2D:  An  Explicit  Two-dimensional  Hydrodynamic  Finite 
Element  Code  with  Interactive  Rezoning  and  Graphical  Display”,  LSTC  Report  1004, 
Livermore  Software  Technology  Corporation,  1990 


UNCLASSIFIED 


INTERNAL  DISTRIBUTION 
DREY  -  TM  -  9822 


1  -  Deputy  Director  General 

1  -  Chief  Scientist 

6  -  Document  Library 

1  -  Head,  Weapons  Effects  Section 

1  -  Dr.  G.  McIntosh  (author) 

1  -  Dr.  D.  Nandlall 

1  -  Dr.  K.  Williams 

1  -  Mr.  M.  Szymczak 

1  -  Mr.  G.  Pageau 


UNCLASSIFIED 


EXTERNAL  DISTRIBUTION 
DREV  -  TM  -  9822 


1  -  DRDCIM 
1  -  DRDCIM  (unbound) 

1  -  DRDB 
1  -  DSAL 
1  -  DSAL- 2 
1  -  DSAL-4 

Attn:  Mr.  David  Saint 
1  -  DAES 
1  -  D  AES-3 
1  -  DFTEM-1 
1  -  DRES 

Attn:  Dr.  C.  Weickert 

1  -  National  Library  of  Canada 

1  -  CISTI 
1  -  DTIC 

1  -  Livermore  Software  Technology  Corporation 

Attention:  Dr.  J.O.  Hallquist 
2876  Waverley  Way 
Livermore,  California,  USA  94550 

1  -  US  Army  Research  Laboratory 

AMSRL-WM-MF 
Attention:  Dr.  A.M.  Rajendran 
Aberdeen  Proving  Ground 
Aberdeen,  Maryland,  USA  21005 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  FORM 
(Highest  classification  of  Title,  Abstract,  Keywords) 


DOCUMENT  CONTROL  DATA 

1 .  ORIGINATOR  (name  and  address)  2.  SECURITY  CLASSIFICATION 

Defence  Research  Establishment  Valcartier  (Including  special  warning  terms  if  applicable) 

2459  Pie-XI  Blvd  North  Unclassified 

Val-Belair,  Qc  G3J  1X5  Canada 


3.  TITLE  (Its  classification  should  be  indicated  by  the  appropriate  abbreviation  (S,  C,  R  or  U) 

The  Johnson-Holmquist  Ceramic  Model  as  used  in  LS-DYNA2D 

4.  AUTHORS  (Last  name,  first  name,  middle  initial.  If  military,  show  rank,  e.g.  Doe,  Maj.  John  E.) 

McIntosh,  G. 

5.  DATE  OF  PUBLICATION  (month  and  year) 

6a.  NO.  OF  PAGES 

6b.  NO.  OF  REFERENCES 

1998 

21 

3 

7.  DESCRIPTIVE  NOTES  (the  category  of  the  document,  e.g.  technical  report,  technical  note  or  memorandum.  Give  the  inclusive  dates  when  a  specific 
reporting  period  is  covered.) 
technical  memorandum 


8.  SPONSORING  ACTIVITY  (name  and  address) 


9a.  PROJECT  OR  GRANT  NO.  (Please  specify  whether  project  or  grant)  9b.  CONTRACT  NO. 
2bb15,  Numerical  Modelling  of  Ballistic  Events 


10a.  ORIGINATOR'S  DOCUMENT  NUMBER  10b.  OTHER  DOCUMENT  NOS 

DREV-TM-  9822 

N/A 

1 1 .  DOCUMENT  AVAILABILITY  (any  limitations  on  further  dissemination  of  the  document,  other  than  those  imposed  by  security  classification) 

E3  Unlimited  distribution 

LJ  Contractors  in  approved  countries  (specify) 

1  1  Canadian  contractors  (with  need-to-know) 

l~l  Government  (with  need-to-know) 

□  Defense  departments 

I  I _ Other  (please  specify) _ _ _ 

12.  DOCUMENT  ANNOUNCEMENT  (any  limitation  to  the  bibliographic  announcement  of  this  document.  This  will  normally  correspond  to  the  Document 
Availability  (11).  However,  where  further  distribution  (beyond  the  audience  specified  in  11)  is  possible,  a  wider  announcement  audience  may  be 
selected.) 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  FORM 
(Highest  classification  of  Title,  Abstract,  Keywords) 

dcd03e 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  FORM 
(Highest  classification  of  Title,  Abstract,  Keywords) 


dcd03e 


